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ABSTRACT 

We investigate the nature of the mass-metallicity (M-Z) relation for long gamma-ray burst (LGRB) 
host galaxies. Recent studies suggest that the M-Z relation for local LGRB host galaxies may be 
systematically offset towards lower metallicities relative to the M-Z relation defined by the general star 
forming galaxy (SDSS) population. The nature of this offset is consistent with suggestions that low 
metallicity environments may be required to produce high mass progenitors, although the detection 
of several GRBs in high-mass, high-metallicity galaxies challenges the notion of a strict metallicity 
cut-off for host galaxies that are capable of producing GRBs. We show that the nature of this reported 
offset may be explained by a recently proposed anti-correlation between the star formation rate (SFR) 
and the metallicity of star forming galaxies. If low metallicity galaxies produce more stars than their 
equally massive, high-metallicity counterparts, then transient events that closely trace the SFR in a 
galaxy would be more likely to be found in these low metallicity, low mass galaxies. Therefore, the 
offset between the GRB and SDSS defined M-Z relations may be the result of the different methods 
used to select their respective galaxy populations, with GRBs being biased towards low metallicity, 
high SFR, galaxies. We predict that such an offset should not be expected of transient events that do 
not closely follow the star formation history of their host galaxies, such as short duration GRBs and 
SN la, but should be evident in core collapse SNe found through upcoming untargeted surveys. 
Subject headings: gamma rays: bursts — galaxies: star formation 



1. INTRODUCTION 

Investigations of the environments in which gamma- 
ray bursts (GRBs) occur has long been an important 
path to understanding the nature of their progenitors, as 
different progenitor models have traditionally predicted 
distinct GRB host galaxy populations. The observations 
of long (LGRB) host galaxies, made possible through X- 
ray localizations of GRB afterglows, show that they tend 
to be bluer, fain ter, and later type than galax ies at 
similar redsh i fts (IFruchter et al.lll999i: iChary et alj|2002t 
Bloom et all 120021: iLe Floc'h et al.l f2003t iTanvir et all 
20041: IFruchter et all 120061: ICastro Ceron et al.l 120061 ). 
They also tend to have higher specific star forma- 
tion than typical star- fo rming galaxies (|Charv et all 
12001 iBerger et al] [200l iChristensen et al.l I2004J) and 
their redshift distribution tends to broadly track the 
overal l cosmic sta r formation rate (SFR) of the uni- 
verse (Bloom 2003; Firmani et al . 2004; Nataraia n et al.l 



20051: iJakobsson et all 120061: iKocevski fc Liangl 120061 : 
Guctta & Pira nl 120071 ). These results have bolstered 
the connection between GRBs and the death of mas- 
sive stars, which is now well-established at low redshift 
(z < 0.3) through the association of GRBs with broad 
lined SN lc events (for a review, see iWooslev fc Blooml 
120061) . 

Spectroscopic observations of host g alaxies have also 
show n that they tend to b e met al poor dProchaska et all 
20041: [SoTlerman et al] 120051: IFruchter et al.l l200e 



Modiaz et al] 120061: [Stanek et all 120061: iThone et all 



20071: iWiersema et al.l (2007: M argutti et all 120071 ). A 
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detailed comparison between the metallicity at the sites 
of broad-lined SN Ic that have been associated with 
GRBs and the site of SN Ic with no detected gamma-ray 
emission found that the chemical abundance of SN-GRB 
hosts were syst ematically lower tha n the hosts of SN 
without GRBs iModiaz et all ((20081 ) . This metallicity 
difference raised the possibility of a upper limit to the 
metallicity of a galaxy that can produce a GRB. 

Recently, lHan et al ((2010) and iLevesque efall (|2010l ) 
compared the mass-metallicity (M-Z) relation for long 
GRB host galaxies to samples from the Sloan Digital Sky 
Survey (SDSS) rep resentative of the gene ral star-forming 
galaxy population dTremonti et al.|l2004D . Using a small 
sample of 5 host galaxies . lHan et al (2010) found that the 
metallicities of the host galaxies tended to fall below the 
low r edshift M-Z rela t ion de fined by SDSS catalog. Like- 
wise, ILevesque et al.l ((2010) compared a much broader 
sample of LGRB host galaxies and found a similar off- 
set, with LGRB host galaxies exhibiting lower metallici- 
ties compared to SDSS galaxies of similar masses. 

The relative nature of this metallicity offset for a 
given host mass, along with a small but growing num- 
ber of hig h-mass, high-m et allicit y host galaxies pre- 
sented by ILevesque et al.l (I2010D and o ther studies 
(iFvnbo et all 12006): iProchaska et ail 120071 : IFvnbo et all 
120081: iChen et al] 12008ft. challenge the notion of a sharp 
meta llicity threshold (jModiaz et all20 06: Koc evski et all 
12009ft for the host galaxies that are capable of producing 
LGRBs. Moreover, without such a strict low metallicity 
criteria for the generation of GRB progenitors, the reason 
of the preference of GRBs to occur in environments with 
relatively lower chemical enrichment remains unclear. 

In this letter we examine the mass-metallicity rela- 
tion for GRB host galaxies by investigating the effects of 
a recently proposed anti-correlation between a galaxy's 
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Fig. 1. — An attenuated power-law model for the SFR-Z relation 
for values of log(M*/M Q ) between 9.0 < log(A/*/M ) < 11.5 at 
intervals of 0.5 log(M* ./ Mq) . The median SFR increases sharply 
as a function of falling metallicity for low mass galaxies, but no 
such relation is observed for their high mass counterparts. 

star formation rate and its metallicity. By examin- 
ing star forming gala xies from a number of catalogs, 
iMannucci et al.l (|2010f ) found a trend of increasing SFR 
as a function of decreasing metallicity in low-mass galax- 
ies. While this trend is prominent in low-mass galaxies, 
the significance of the anti-correlation decreases sharply 
with increasing galaxy mass, almost disappearing com- 
pletely for high-mass galaxies. 

We show that such a connection between the SFR 
and chemical enrichment of a galaxy w ould naturally ex - 
plai n the systematic offsets reported by lHan et all (|2010l ) 
and iLevesque et al.l (|2010f ). If low metallicity galaxies 
produce more stars than their equally massive, high- 
metallicity counterparts, then it would be expected that 
transient events that closely trace the star formation his- 
tory of their host galaxies would be more likely to occur 
in these high SFR, low-metallicity galaxies. We test this 
hypothesis by modeling the M-Z relation for LGRBs us- 
ing a combination of data from the SDSS-DR7 catalog 
and model prescription for the SFR-Z relation in Section 
2. We compare our model predictions to published host 
mass data in Section 3 and we discuss the implications 
of our results in Section 4. 

2. MODEL PRESCRIPTIONS 

In order to investigate the nature of the metallically 
offset in the M-Z relation for LGRB host galaxies, we 
must consider the difference between the probability dis- 
tributions for the metallicity of a galaxy drawn from 
the general star forming population and a galaxy that 
is likely to harbor a LGRB. For the general star form- 
ing population, the metallicity probability distribution 
P(Z) for a galaxy of a given stellar mass M* is simply 
the normalized metallicity distribution <f>(Z, M*) of all 
star forming galaxies at that mass. Therefore, the most 
likely metallicity of a galaxy drawn randomly from the 
star forming population will reflect the peak of <fi(Z, M*). 
The most likely metallicity of a galaxy selected because 
of the occurrence of a LGRB is quite different. For 
host galaxies, any metallicity dependance on the prob- 
ability of the transient event, Pqkb(Z, M*), must be 
taken into consideration. In such a case, the resulting 
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Fig. 2. — The normalized metallicity distribution (f>(Z) (blue 
histogram) for the general star forming population and the SFR 
weighted metallicity distribution Pn OB t(Z, M*) distribution (red 
histogram) between 9.9 < log(A/*/Af Q ) < 10.1. The peak of the 
-Pliost^) distribution is distinctly offset to lower metallicities be- 
cause of the greater weight given to low metallicity galaxies due to 
their higher SFR compared to their higher metallicity counterparts. 

LGRB metallicity distribution would reflect the product 
of the metallicity distribution for all star forming galaxies 
and this metallicity dependent probability distribution, 
iW (Z) ~ 4>{Z, M*) * P LGRB {Z, M*). 

The metallicity dependance on the probability of the 
transient event can represent a number of different the- 
oretical considerations. It could reflect, for example, 
the metallicity dependance on the likelihood that a star 
would possess the physical conditions required to pro- 
duce a GRB. Conversely, it could even reflect an ob- 
servational bias such as the hypothetical probability of 
optically detecting the event in galaxies with metallicity 
dependent degrees of dust obscuration. In this letter we 
focus on the metallicity dependance of a galaxy's star for- 
mation rate as the primary contributor to Plgrb (Z, M+). 

Because the production (and death) of massive stars is 
closely linked to the ongoing star formation in a galaxy, 
we assume that a higher SFR in a galaxy will result in a 
greater likelihood for the occurrence of a LGRB. There- 
fore we set Plgrb(^, M*) ~ SFR(Z, M+). Here we are 
explicitly assuming that an equal fraction of stars being 
produced, in both low and high SFR galaxies, contribute 
to the production of massive LGRB progenitors. 

In order to model the M-Z relation for GRB host 
galaxies, we must convolve the <fi(Z) and Plgrb (Z) prob- 
ability distributions and find the peak of the result- 
ing PHost(-Z) distribution for a range of galaxy masses. 
This requires adopting prescriptions for both <p(Z) and 
SFR(Z) as a function of galaxy mass. We can mea- 
sured cb(Z) directly through the use of data from the 
MPA/JHU catalog from SDSS-DR70. This catalog in- 
cludes stellar masses and emission line estimates for over 
927552 galaxies, providing a wealth of information re- 
garding galaxy demographics in the local Universe. The 
metallicity values presented in the catalog reflect oxygen 
abundan ce estimates based on a statistical method dis- 
cussed in lTremonti et all uj2004T ) . We also utilize the cat- 

3 http:/ /www. mpa-garching.mpg.de/SDSS 

4 http://www.sdss.org/dr7 
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Fig. 3. — The M— Z relation for the general star forming population between 0. 005 < z < 0.30 , with the color coding representing the 
absolute SFR. The blue diamonds represent the median M— Z relation presented in Tremonti et al. (2004) with the dashed and dotted lines 
representing the contours which include 68% and 95% of the data, respectively. The red circles represent our predicted median GRB host 
galaxy metallicity as a function of log(M*), effectively producing an SFR weighted M-Z relation which is shifted to lower metallicitics. 
The vertical error bars represent the ler spread in the unde rlying -Pjiost (Z) distribution. The blue circles represent the mass and metallicity 
of GRB host galaxies selected from Lcvcsauc et al. (20f0) for which a conversion between the KK04 and T04 metallicity diagnostics was 
possible. 



alog's stellar masses, obtained through broad-band fits 
SDSS photometry and SFR estimates based on the tech- 
nique discussed in iBrinchmann et"aL1 (|2004fl . 

For the purposes of our model, we selected emission 
line galaxies with a redshift cut of 0.005 < z < 0.30 
for which stellar mass estimates were possible. For the 
remaining galaxies, we set a threshold to the signal-to- 
noise ratio (S/N) of Ha of S/N > 5 in order to insure 
reliable metallicity estimates. Metallicities for AGNs are 
not reported in the catalog and hence are not included in 
our sample. Our final sample contains 155564 galaxies. 

For an analy t ic des cription of SFR(Z), we turn to 
iMannucci et al.l (|2010l ). where the authors quantified a 
trend between SFR and metallicity using the SDSS-DR7 
MPA/JHU catalog at redshifts between 0.07 and 0.30. 
They find that for all log(M*) < 10.7, SFR increases with 
decreasin g metallicity at constant mass (e.g. see their 
Figure 1. IMannucci et al.l (|2010P i). They provide a 4th 
order polynomial fit to the observed SFR-Z relation, al- 
though for the purposes of our analysis we have reworked 
their expression into an attenuated power-law in order to 
extrapolate their original fit below log(SFR) < — 1 and 
above log(SFR) > 1. Our analytic model for SFR(Z) 
is presented in Figure Q] for several different values of 
log(M*/M ) between 9.25 < log(M*/M ) s; 11.35. 
We adopt this empirical model as the prescription for 



SFR(Z) as a function of M*/M . 

3. RESULTS 

Using the subset of data from MPA/JHU catalog, 
we can obtain <j>(Z) for effectively a constant mass by 
selecting all galaxies within a narrow M+/Mq range. 
For example, in Figure [2] we show the normalized cj>(Z) 
(blue histogram) selected for all galaxies between 9.9 < 
log(M*/M Q ) < 10.1 cf>(Z), totaling 20450 galaxies. The 
peak of this distribution represents the median metalli- 
cally for a star forming galaxy in the MPA/JHU catalog 
for log(M*/M ) ~ 10.0 ±0.1. We can then calculate the 
median metallicity of a galaxy likely to host a GRB by 
weighting (j>(Z, M*) by SFR(Z,M*) and finding the peak 
of the resulting probability distribution Pn os t(Z, M+). 
The normalized Pn os t(Z, M*) distribution is shown in 
Figured] (red histogram). The peak of the PH OS t(Z, M*) 
distribution has shifted to lower metallicities because of 
the greater weight given to low metallicity galaxies due 
to their higher SFR compared to their higher metallicity 
counterparts. 

We calculate the median of the 4>{Z) and Pnost(-^) 
distributions for a range of log(M^/M ) values between 
8.0 < log (M+/M®) < 11.0 and plot the results in Figure 
[H The red circles represent the predicted median GRB 
host galaxy metallicity as a function of log(M*/M ), ef- 
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Fig. 4. — The SFR vs. M^/Mp, fo r all 155564 galaxies in our SDSS sample along with 6 galaxies presented in Levesque et al. (2010) and 
2 galaxies from (Savaglio et al. 2009) that fell below z < 0.3. The color gradient represents metallicity and the expected trend of increasing 
SFR and metallicity as a function of galaxy mass is clearly evident. All 8 of the low redshift LGRB hosts under consideration lie to the 
upper left of the SFR- M+/M Q trend. 



fectively producing an SFR weighted M-Z relation that 
is shifted to lower metallicities compared to the M-Z 
relation for the general star forming galaxy population. 
The vertical error bars represent the la spread in the 
underlying Pn os t(Z) distribution. The offset between the 
SFR weighted and general M-Z relations is roughly ~ 0.3 
at log(M + /M ) - 9.0 and - 0.15 at log(M*/M ) - 11.0. 
The decrease in the offset between the SFR weighted 
and general M-Z relations is con sistent with the trend 
reported bv lLevesaue et all (|2010f ) in which the metallic- 
ities of high-mass host galaxies are more consistent with 
the general M-Z relation compared to their low-mass 
counterparts. 

We compare our modeled SFR weighted M-Z relation 
to the mass and met allicity of a sub s et of GRB host 
galaxies selected from iLevesque et al.l (|2010[ ). Because 
ILevesque et alJ (|2010f ) report their metallicity values in 
the KK04 system, we convert their reported value to 
the T04 system using t he co nversion tables provided in 
IKobulnickv fc Kewlevl (|200l . Since this conversion is 
not defined for all metallicities, we selected a subsam- 
ple of 4 host galaxies below z < 0.3 that fell within the 
metallicity range for which this conversion was possible. 

Although the predicted SFR weighted M-Z relation is 
more consistent with the data than to the general M-Z 
when considering the errors in both mass and metalli- 
cally, all four galaxies still nonetheless systematically fall 
below our model. We return to the implications of this 



discrepancy in the next section. 

4. DISCUSSION 

The nature of the low metallic ity offset in the GRB de- 
fined M-Z relation reported bv ILevesque et ail (|2010f ) is 
rather intriguing. The observed shift supports previous 
spectroscopic evidence that GRBs tend to preferentially 
occur in low metallicity environments, although the rela- 
tive nature of the offset, in additio n to several high metal - 
licit y host galaxies repo rted by ILevesque et al.l (|2010D 
and (| Savaglio et al. 2009) , suggests that host galaxies do 
not necessarily adhere to a strict metallicity "cut-off. 

The possibility of a correlation between the chemical 
enrichment and the rate of star formation in a galaxy may 
begin to explain the nature of this relative offset, with- 
out the need of such a cut-off in terms of absolute metal- 
licity. If low metallicity galaxies tend to produce more 
stars than their high metallicity counterparts for a given 
stellar mass, then the likelihood of detecting a transient 
event linked to a massive, short lived progenitor from 
within a low metallicity galaxy will increase accordingly. 
Therefore the nature of the offset of the GRB defined 
M-Z relation towards lower metallicity could largely be 
explained as a bias towards finding transient events in 
galaxies with relatively higher SFRs. 

Several studies have already indicated that LGRBs do 
tend to occur in galaxies with higher specific star forma- 
tion than typical star forming galaxies at similar redshifts 
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(iCharv et al.ll2002t iBerger et al.|[200llChristensen et al.l 
12004 ISavaglio et al.ll2009D . We illustrated this in Figure 
HI where we plot the SFR vs. M*/M for all 155564 
galaxies in our SPSS sam p le alo ng with 6 galaxies pre- 
sented in ILevesque et al.l (|2010l ) and 2 galaxies from 
(jSavaglio et al.ll2009D that fell below z < 0.3. The color 
gradient represents the metallicity of the SDSS galaxies 
and the expected trend of increasing SFR and metallicity 
as a function of galaxy mass is clearly evident. All 8 of 
the low redshift LGRB hosts under consideration occupy 
the upper left of the SFR- Z-M+/M Q plot, having SFRs 
that are substantially larger than the mean SDSS sam- 
ple for their galaxy mass. These few galaxies illustrate 
the degree to which these host galaxies have SFRs that 
exceed the star formation activity of similar galaxies. 

A bias towards finding transient events in high SFR 
galaxies alone does not explain why LGRBs occur in 
small to intermediate size galaxies as opposed to the most 
massive galaxies, which typically have much higher ab- 
solute rates of star formation. To understand the mass 
distribution of LGRB host galaxies, we must consider 
the number of galaxies as a function of mass (the galaxy 
mass function) . Although high mass galaxies produce far 
more stars than their low mass counterparts at all red- 
shifts, the number of low mass galaxies in the Universe 
exceeds that of high mass galaxies by several orders of 
magnitude. As a result, the SFR weighted galaxy mass 
function, i.e. the galaxy mass at which most of the star 
formation is occurring at a given redshift, peaks at inter- 
mediate masses. Therefore the fact that we do not see 
LGRBs occurring predominately in the most massive, 
highly star forming, galaxies is not in contradiction to 
the suggestion that transient selected galaxies are biased 
tow ards high SFR environ ments . 

In lKocevski et all ()2009l ) we compared the mass distri- 
bution of GRB host galaxies to the SFR weighted galaxy 
mass function. Using a similar argument as we have in 
this letter, we modeled the SFR weighted galaxy mass 
function for a variety of redshifts by convolving prescrip- 
tions for the the galaxy mass function <P{M^/Mq) at 
a given redshift with the SFR as a function of stellar 
mass. We found that the SFR weighted mass function 
at z = 1 peaks at roughly log(M*/MQ) ~ 10.3. The 
measured median stellar mass o f LGRB host galaxies re- 
ported by Savagli o et al.l (|2009l ) for the same redshift is 
slightly lower, roughly log(M*/M Q ) ~ 9.3. 

We interpreted this discrepancy as evidence for a 
metallicity cut-off at which the SFR weighted galaxy 
mass function would be truncated due to the relationship 
between a galaxy's mass and metallicity, resulting in a 
host mass distribution which would be shifted to lower 
masses. Such a metallicity cut-off may not be necessary 
if in fact low mass, low metallicity galaxies produce more 
stars than their high metallicity counterparts. In such a 
scenario, the SFR weighted galaxy mass function would 
be shifted to lower masses without the need of a metal- 
licity cut-off due to the additional weight given to low 
mass, low metallicity galaxies. The effects of a SFR-Z 
relation for low mass galaxies was not taken i n to acc ount 
in the analysis performed by iKocevski et al.l ( 2009) and 
a full examination of the effects of the correlation on the 
resulting joint mass and metallicity distributions will be 
reserved for a future paper. 

Although we have shown that a metallicity dependance 



on the SFR of low mass galaxies can in fact produce 
an M-Z relation offset to lower metallicities, th e mass 
and metallicity data taken from ILevesque et al.l (j2010f ) 
nonetheless systematically fall below our modeled M-Z 
relation. Within the context of our model, this discrep- 
ancy may be evidence for a further bias in the probability 
of LGRBs to occur in relatively low metallicity environ- 
ments beyond the metallicity dependent SFR that we 
have considered. A sharp decrease in Plgrb(Z) with 
increasing metallicity due to a metallicity dependance 
on the physical conditions required to produce a GRB, 
would be one such effect. In the framework of the collap- 
sar model, low metallicity progenitors would retain more 
of thei r mass due to the reduc t ion of line-driven stella r 
winds (IKudritzki fc Pulsl IIM IVink fc de Koterl [2005T) . 
and hence preserve more of their a ngular momentum and 
stellar mass at the time of collapse ()Yoon fc Langer|[2005l : 
lWooslev~fc Hcgcr 2006) . This mechanism may be crucial 
to the production of collimated emission that is thought 
to be required to produce LGRBs. 

Therefore, although the observations of several high 
mass, high metallicity LGRB host galaxies may rule out 
a strict metallicity "cut-off' , t he discrepancy betwee n the 
host galaxy data measured by ILevesque et all (|2010l ) and 
our modeled M-Z relation may point to further metallic- 
ity biases in the LGRB sample. These biases need only 
to reflect a decreasing likelihood of a GRB occurrence as 
a function of increasing metallicity in order to provide a 
mechanism to further shift the LGRB defined M-Z rela- 
tion to lower metallicities. Such biases would compound 
the effect introduced through the simple metallicity de- 
pendent SFR prescription that we have assumed in our 
model and as such there would be no need for a strict 
metallicity "cut-off" to explain the relative preference for 
LGRBs to occur in low metallicity environments. 

Finally, our work predicts that core collapse SNe found 
through untargeted surveys such as the Palomar Tran- 
sient Factory (PTF), Pan-Starrs, and the Large Synoptic 
Survey Telescope (LSST) should also define an M-Z re- 
lationship that is shifted to lower metallicities compared 
to the general star-forming galaxy population. Such an 
effect should not be expected of transient events that do 
not closely follow the star formation history of their host 
galaxies, such as short duration GRBs and SN la events. 
If subsequent observations of SN II and Ibc host galaxies 
found through these surveys do not show such an offset, 
then the shift in the GRB defined M-Z relation would 
necessarily be a physical effect intrinsic to GRBs and not 
due to the SFR-Z selection bias that we propose. 

5. CONCLUSIONS 

A growing body of spectroscopic evidence has made 
it clear that GRBs tend to preferentially occur in host 
galaxies that exhibit lower metallicities compared to the 
general star-forming galaxy population. Recent stellar 
mass and metallicity estimates made by ILevesque et al.l 
(2010) indicate that host galaxies follow a M-Z rela- 
tion that is offset to lower metallicities c ompared to the 
genera l star-forming galaxy population. ILevesque et al.l 
(2010) find that this offset, while pronounced for low 
mass galaxies, becomes negligible for two GRBs found 
to occur in high-mass, high-metallicity host galaxies. 

We show that the nature of this reported offset may 
be explained by a recently proposed anti-correlation be- 
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tween the SFR and the metallicity of star forming galax- 
ies. If low metallicity galaxies produce more stars than 
their equally massive, high-metallicity counterparts, then 
it would be expected that transient events that closely 
trace the star formation in a galaxy would be more likely 
to be found in these low metallicity, low mass galax- 
ies. We quantify this hypothesis by modeling the M- 
Z relation for GRBs. We calculate the SFR weighted 
metallicity distribution for star forming galaxies by us- 
ing metallicity data from the MPA/JHU catalog as part 
of SDSS-DR7, as well as a model prescription for the 
relationsh ip between a galaxy's SFR and metallicity re- 
ported bv lMannucci et all ()2010D . We find that the SFR 
weighted metallicity distribution peaks at significantly 
lower metallicity compared to the metallicity distribu- 
tion for the general star forming galaxy population at a 
given stellar mass. 

Therefore, the offset between the GRB M-Z relation 
and the properties of the general star forming population 



could simply be due to the methods used to select their 
respective galaxy populations, with GRBs being biased 
towards low metallicity, high SFR, galaxies. Although 
our model provides a framework to explain the nature 
of this low metallicity offset, we find that further biases 
beyond a metallicity dependent star formation may be 
required to account for the full degree of the discrepancy 
between the GRB and SDSS defined M-Z relation. 
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